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We theoretically investigate the response of a two component Fermi gas to vector potentials 
which couple separately to the two spin components. Such vector potentials may be implemented in 
ultracold atomic gases using optically dressed states. Our study indicates that light-induced gauge 
potentials may be used to probe the properies of the interacting ultracold Fermi gas, providing, 
amongst other things, ways to measure the superfluid density and the strength of pairing. 

PACS numbers: 67.85. Lm,47.37.+q 



I. INTRODUCTION 

Experimental progress in the trapping and manipula- 
tion of ultracold Fermi gases has led to new regimes of 
study of superfluid two component Fermi systems. This 
is thanks to the unprecedented control over the micro- 
scopic properties of many-body systems which ultracold 
Fermi gases offer. Examples of this microscopic control 
include the modification of the interaction strength via 
Feshbach resonances [11, Q , the tuning of the density im- 
balance of spin up and spin down particles 0, |3] and 
the ability to impose controllable lattice potentials us- 
ing optical lattices In particular the tunability of 
the interaction strength has allowed interesting investiga- 
tions of the continuous crossover from a Bardeen-Cooper- 
Schrieffer (BCS) type superfluid, analogous to the BCS 
state in a superconductor, to a Bosc Einstein condensate 
(BEC) of molecules consisting of bound spin up and spin 
down pairs 

Experiments on ultracold atomic Fermi gases allow 
studies of a variety of physical properties of a fundamen- 
tal interest. Past work has studied the pairing gap as a 
function of the temperature and the interaction strength 
pjj . the transition from superfluid to normal behaviour 
by studying collective modes and the formation of 
vortices in rotating ultracold Fermi gases [T3| . 

Light-induced gauge potentials [l4| offer the possibility 
to study new aspects of the response of superfluid Fermi 
system. This is the issue that we explore in this paper. 
As we describe below, light-induced vector potentials of- 
fer the opportunity to study diverse properties of ultra- 
cold Fermi gases, ranging from the superfluid density to 
properties of the pairs. These different quantities become 
experimentally accessible by varying the time depen- 
dence of the light induced vector potential and by con- 
sidering either the response to a vector potential which 
couples equally (spin-symmetric) or in the opposite way 
(spin-asymmetric) to the two spin species. Specifically, 
we decompose the vector potentials A-|-^ for the two spin 
species into a symmetric component A = ^{A^+Ai) and 
an asymmetric component A A = A^ — A^. We study 



the response of a two-component Fermi gas, working in 
the BCS limit, separately to the spin-symmetric A and 
spin-asymmetric AA vector potentials. 

The paper is organised as follows. In Sec. |ll] we de- 
scribe how the spin-symmetric vector potential A leads 
to a superfluid response. In Sec. |lll] we describe the re- 
sponse to a spin-asymmetric vector potential AA in the 
zero temperature limit, exploring the frequency depen- 
dence in detail. The effects of non-zero temperature are 
described in Sec. lIVI In Sec.|V]we explain how the effects 
we predict can be studied in experiment. Finally, Sec. lVIl 
summarizes our main results. 



II. EFFECT OF A SPIN-SYMMETRIC VECTOR 
POTENTIAL 

As discussed in [l^ [l^ for the bosonic system, the 
response of a one component Bose gas to a vector poten- 
tial is determined by the superfluid density of the Bose 
gas. The light-induced vector potential leads to a kinetic 
energy of the form 2^ {p — A)'^ , which has its minimum 
shifted to p = A. The normal fluid component will seek 
the new minimum of the dispersion and thus come to 
rest at the new equilibrium. The superfluid component 
on the other hand is unaffected by the shifting of the 
minimum in the dispersion and will continue its initial 
state which now no longer is the equilibrium state. The 
phenomenology of a two component Fermi gas subjected 
to a spin-symmetric vector potential is the same as that 
of a one component Bose gas subjected to a vector poten- 
tial. For then the two components of the superfluid are 
affected equally by the vector potential and for the entire 
system there is a new steady state. As in the Bose case, 
only the normal part of the superfluid will come to rest 
in the new steady state. The superfluid component will 
not relax to this new equilibrium state, but continue in 
its metastable state, thereby allowing one to distinguish 
between normal and superfluid densities. 

As in the Bose case, the low-frequency response to A is 
simply determined by the superfluid density of the Fermi 
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gas, with a mean fermion number current density 



and the gap is given by [ij 



j^-ij, + j,) = -—A 



(1) 



where ps is the superfluid density, defined such that 
its maximum value is equal to the number density of 
fermions. We will therefore concentrate on the response 
to AA. As we shall show, in some regimes the response 
to AA vanishes, under which conditions the response can 
be determined from the response to A alone and is a mea- 
sure of the superfluid density p^. 



III. EFFECT OF A SPIN-ASYMMETRIC 
VECTOR POTENTIAL 

A. Static vector potential 

We consider the BCS Hamiltonian for a system of spin 
i fermions with vector potentials A^ and A^ coupling to 
the spin up and down particles respectively. It is given 

by 
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+yc|(r)c|(r)c^(r)cf(r) 



(2) 



Here Ca\r) are fermionic annihilation (creation) oper- 
ators for spin a and V is the strength of the contact 
interaction V{r — r') = VS{r — r'). V is related to the 



s-wave scattering length as via V 



n 



We will 



be considering attractive interactions, so < 0. For a 
uniform vector potential with A^ = —A^_ ~ A the physi- 
cal properties of the system are the same as for A = 0. In 
order to show this, perform the following transformation 
on the creation and annihilation operators. 



(3) 



Here and in the remainder of this paper we set h = 1. 
Upon inserting equations ([3]) into equation ^ we recover 
the same Hamiltonian but now for the operators 
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(4) 



For A 



t 



= A the mean- field ground state quan- 
tities pAix) and Aa{x) are the same as the quantities p 
and A for A^- = 0. The densities are given by 



(5) 



AA{r) = Feff (ci(r)ct(r)) = (d^(r)dt(r)e'(^^+^^'-'- 
=A(r) = Ao (6) 



since A is homogeneous. We have introduced an ef- 
fective potential Vi^fr to regularise the divergent nature 
of the contact interaction [13] • This shows that for 
Af = ~Ai — A the ground state properties are indeed 
independent of A. 

This is to be contrasted with the result for Af = Aj^. 
In that case the ground state of the system is a state in 
which the phase of the gap is spatially varying. This can 



be seen from equation 
into equation ([6]) we obtain 



If we insert A 



t 



AA{r) = K=ff (ci(r)ct(r)) = V,ff (^d^{r)d^{ry^^^+^^ 

= Aoe2'^-. (7) 



showing that the phase of the gap is (p = 2\A-r and there- 
fore spatially varying for nonzero A. Although the gap is 
spatially varying in the groundstate, the gauge-invariant 
supercurrent density remains zero. The density also re- 
mains unchanged, as can be seen from equation ([5]). 



B. Time-dependent vector potential 



Given that the application of a static vector potential 
A with A^ = —A^ = A does not have an effect on the 
ground state of the system we consider a time-dependent 
vector potential A. In the analogous electromagnetic sys- 
tem 7^ corresponds to the presence of an electric 



field [18[. However, since we are considering A-f = — Aj, 
this would correspond to a Fermi system where the spin 
up and spin down particles have opposite charges. An 
example would be a spin polarised electron-positron su- 
perfluid where the spin up particles are positrons and the 
spin down particles are electrons. 

We now treat the application of time-dependent vec- 
tor potentials Aa{t) within linear response. Write the 
Hamiltonian as 



H = Ha + 5H{t). 



(8) 



Hq is now the mean-field BCS Hamiltonian in the ab- 
sence of a vector potential which we consider in second 
quantised form in momentum space. It is given by 
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h.c. 



(9) 
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A is determined self-consistcntly by the equation [T] 

T7 = ^Eb:--) (10) 



V 2Ld^\Ek 

k ^ 



where L is the size of the system. To linear order in A, 
SH is given by 



4.iWcfc4W ) • (11) 



The response of an observable M is then given by the 
Kubo formula IH 



d{M{t)) = J f dr 



M{t),SH{T) 



(12) 



We consider Aa{t) = Ao-e'"* as the time-dependent per- 
turbation, and study the current response. As we are 
applying a spin-asymmetric vector potential we look at 
the spin current AJ = J-|- — Jj,. Linear response of the 
current involves two contributions. First, there is the 



contribution from the Kubo formula, equation (|12p , with 
M replaced by the spin current operator in the absence 
of a vector potential A 

^'^(*) = E i i^Ui^Mt) - cliit)ck.l{t)) . (13) 



Secondly, there is a contribution arising from the redef- 
inition of the spin current for nonzero A. Applying the 
vector potential A changes the current operator to the 
gauge invariant current operator 



fc A, 



(14) 



This means that definition of the current changes such 
that 



J ->J -V— 

Ja -Ja 2^ C^^c 



Ck,, 



(15) 



Inserting equations (jTT]) , ((T^ and ([T5|) into equation 
we obtain for the spin current in linear response 
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(AJ) {t) = -J2- 
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E-(Ct(^)^fc^t(i)-44(^)^fc4(i)) , 

k 

) ■ (16) 



The term e"''^*""^^ serves to regularise the expression 
and physically means that the perturbation is gradually 
switched on starting at T = — oo. We will eventually take 
the limit 77 — ?> 0. 



To evaluate the expressions in equation (jl6p we expand 
the operators Ck^a in terms of the quasiparticle operators 
dfc o- via [11 



where Uk 



Ck,t 
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Uk Vk 
Vk -Uk 



and Vk 



-T^ with efe = 1^ - ^ 



and iJfe = y^e| + A2. The time dependence of the oper- 
ators ak,a is given by al_^{t) = a^ ^e'^''* and ak^c,{t) = 



(^k,a& Inserting these gives the response 

(^) (^) -A f J2 -ke-^'-'-^^ulvl (2/(i?,) - 1) 
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(k ■ A^{t) - k ■ A^ir)) 
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f{Ek)^]-{A^{t)-A^{t)), 



Ek J m 



(18) 



where / is the Fermi function. Taking A-^ — A^ to point 
along the z-axis and converting the sum to an integral 

(17) via X]fc ^ (^) / '^'^^ ^^^"^ Fourier transforming gives 
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for d = 3 
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1 1 



2Ek + uj + b] 2Ek — uj — irj 



1 / L 
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TO \ 27r 

(^t(^) -^i(^)) 
= x(w)AA(w) 



(19) 
(20) 



where x{uj) is the response function for AJ along the 
z-direction. In the zero temperature hmit x(w) can be 
expressed as 



xH-^(^) ^(2to)^/^A3/^ 



dx- 
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(21) 



find an analytic expression for x('^)- First expand the 
integral in equation (PT|) for w ^ 1 



dx- 
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(.T2-a2)2 + l \ 2^(a;2-a2)2 + i+cD + i7? 
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X^ 



a2)2 + 1)3/2 



(^2 



dx 



((a;2 - a2)2 + 1)5/2 



+ 0(^4). (22) 



Now note that for T = the a;° part of this integral can- 
cels the second term in eq. (|2ip coming from the redefi- 
nition of the current. The integral in the second term in 
eq. (EU can be obtained from [H ^^i_^_^i^2j^1^i^3/2 = 
K{d) 2(c^'^-b'^) ^ E{d) by parametric differentiation 

with respect to c. Here d is given by d = ^yi^^^ ^^"^ 
K{x) and E{x) are complete elliptic integrals of the first 
and second kind and are defined as 



Defining r 



K{x) = 
E{x)^- 



1 



x'^ sm 



1 — a;2 sin 



(23) 
(24) 



%/2(a* + i; 



Y/-r~ we finally obtain the follow- 



ing for the zero temperature response function in the 



limit 3 < 1 



where = and cj = In the limit (2; ^ 1 we can 



^^'^>'~m^\2^) 3^^ 24(a4 + 1)3/4 



i/i (r) ((a-* + 1)3/2 _^ 3^4(^4 ^ 1)1/2 _ 4^2^^4 _^ ^ 4i;(r)(a4 + 1)0^ 



(25) 



As the BCS approximation we are working in is only valid 
up to order ^ we expand equation (|25p to lowest order 
in - and obtain 

a 

Equation (PS)) shows that in the limit cj — > 0, 
x(w,r = 0) vanishes quadratically with lo. We thus find 
that the low-frequency response is consistent with the re- 
sult in section lTlI Al namely that a static vector potential 



acting in the opposite way on the two spin species has no 
effect on the ground state of the system. 



C. Interpretation of the low- frequency response 

The low-frequency response of the system can be un- 
derstood in terms of currents arising from the time- 
varying polarisation of the Cooper pairs in the Fermi 
gas. A time dependent vector potential A gives rise to 
an effective electric field E via E = Since we are 
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considering = —A^ this gives rise to an electric field 
which couples differently to the spin up and spin down 
particles, such that we obtain 



AE = -{A^ - Ai) = -iujAA. 



(27) 



In the absence of such an effective electric field the s-wave 
Cooper pairs which make up the superfluid are spheri- 
cally symmetric. Application of an effective electric field 
leads to a displacement of the average positions of the 
spin up and spin down particles. This leads to a polari- 
sation of the Cooper pairs just like the application of an 
electric field to a hydrogen atom induces a dipole mo- 
ment and hence a polarisation in the hydrogen atom. A 
time-dependent polarisation results in a current which is 
given by [ill 



J = 



dP 



(28) 



Combining equation ((27)) and equation p8| we obtain 
the polarisation 



(29) 



which implies that is the polarisability 7. According 
to equation (|26p 7 is constant for small lo and propor- 
tional to A~^. The polarisability describes how easily 
the particles in the spherically symmetric s-wave Cooper 
pairs can be displaced with respect to each other. A 
larger gap A means that the Cooper pairs are more 
tightly bound and the system has a reduced polarisabil- 
ity. 

Note that the polarisability we describe here is very 
different from the one described in Ref. H^. In this pa- 
per the polarisability we are considering is the dynamic 
polarisability arising from the individual Coop er p airs in 
the absence of pair breaking effects. In Ref. [2J on the 
other hand the static polarisation of an atomic Fermi 
gas cloud is considered which arises when the two spin 
species are subjected to different potentials. The polari- 
sation predicted in Ref. [2^ occurs only once Cooper pairs 
have been broken. 



D. Beyond the low- frequency limit 

The above analysis was done for -j- <C 1. For larger 
values of w we need to solve equation (|2ip numerically. 

As can be seen from equation (pi]) . for lo < 2 A the re- 
sponse is purely real in the limit 77—5-0. Since the imag- 
inary part of the response function describes the energy 
absorbed by the system [23| this means that no energy 
is absorbed. This is what is expected since the pertur- 
bation does not have enough energy to break pairs. We 
have plotted the response function divided by in the 
range < w < 2A in figure [T] for a particular set of pa- 
rameters. This shows that for a reasonably wide range of 




to/ A 



FIG. 1: (Colour online) Response function x divided by cj^ 
and normalised to the number of particles as a function of ^ 
in the region < oj < 2A. The interaction strength is given 
by kpUs = —1.1 where as — is the s-wave scattering 

length. The parameter a = is given by a = 3.7 and r — 1.0. 



Solid Line: 
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as obtained by numerically integrating the 



expression in equation (|21|l . Dashed line: Expression obtained 
from equation (|25p 



Lo the response is indeed well described by equation (|25p . 
Significant departures occur only for > 0.5. In figurc[2] 
the real and imaginary parts of x(a;) for < a; < 3A 
are plotted. As w — 2A the real part of x diverges as 
(2 A — w)"^/^. At w = 2 A the real and imaginary parts 
of X diverge with as 77"^/^. The energy absorption 
of the system diverges since the perturbation can reso- 
nantly couple to pair breaking processes. For w > 2A 
the perturbation has enough energy to break pairs but 
is no longer on resonance. The energy absorption in this 
regime is proportional to (w^ — 4A^)~^/^. 

The real part of the response for u > 2 A is dominated 
by the contribution from the redefinition of the current, 
see equation ([15]). This is because the perturbation has 
enough energy to break pairs and once pairs are broken 
they behave like non-interacting particles, which, in the 
absence of collisions, will not equilibrate to the new low- 
est energy state. This leads to response 



AJ 



N 
'2m' 



-AA 



(30) 



where N is the total number of atoms. This behaviour of 
non-interacting particles can be seen from equation (jl9p 
by setting A = 0. 



IV. EFFECTS OF NON-ZERO TEMPERATURE 
A. Collisionless regime 

For non-zero temperature there is a current response 
even at zero frequency. This can be understood by not- 
ing that at non-zero temperature there are thermally ex- 
cited quasiparticles and within the approximation we are 
working in thermally excited quasiparticles do not inter- 
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FIG. 2: (Colour online) Real and imaginary part of the re- 
sponse function for T = as a function of for the same 
system as in fig. [T] In the limit ry — >■ there is an infinitely 
sharp peak at oi = 2A. Above u — 2A the response is inde- 
pendent of oj and is given by equation psp . 



act. The important effects of quasiparticle collisions will 
be treated in section llVBI For non- interacting quasipar- 
ticles the only contribution to the current comes from the 
redefinition of the current. Within this approximation, 
the response function then measures the density of un- 
paired quasiparticles, reaching its maximum magnitude 
at T > Tc when all particles are unpaired. For tem- 
peratures ^ ^ 1 where Tc is the critical temperature 
x(a; — ?> 0) ^ exp (— -y-). On the other hand, close to Tc 

we find that xi^ 0,T) - x(w ^ 0,Tc) - \T - T^f 
where within mean-field theory /3, the critical exponent 
for A, is given by f3 — 1/2 The response as a func- 
tion of temperature at = is shown in figure [3l Above 
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FIG. 3: (Colour online) Response as a function of t^I both 
below and above Tc for the same system as in fig. [T] Normali- 
sation is with respect to the response of a noninteracting gas, 
N/2m. The kink in the response function occurs at T = Tc. 



Tc we have A = and hence the whole gas behaves like a 
gas of non-interacting particles whose response is simply 
given by the redefinition of the current, equation (1151) . 
for the full density leading to (|30p . 



B. Quasiparticle collisions 

In the above derivation of the response we have ne- 
glected quasiparticle collisions. At non-zero temperature 
thermally excited quasiparticles arc present in the Fermi 
gas. A thermally excited quasiparticle with momentum 
Pi and energy ei — /i can undergo a scattering process 



Pl+P2-^ Pi 



-P'2 



(31) 



where the other particle with momentum p2 can be sit- 
uated above or below the Fermi surface. This leads to a 
non-zero collisional relaxation rate F. As the tempera- 
ture tends to zero the lifetime of an excited quasiparticle 
diverges due to Pauli blocking. Unless w ^ F we cannot 
neglect quasiparticle scattering. 

In order to obtain a rough estimate for the collisional 
relaxation rate F, we consider the quasiparticle scattering 



rate in the Fermi liquid picture |25| . This is an appro- 
priate estimate for T ^ Tc <C Tp when a sizeable density 
of unpaired quasiparticles is thermally excited. Within 
perturbation theory the quasiparticle scattering rate for 
T <C Tf for a thermally excited quasiparticle with energy 
e — /i = fc^r and for a S function interaction potential 
V = VS^{r — r') cf. equation can be derived along 
the lines of Ref. M- We find 



r = nikpasf ( ^ 



2 2 



1 



(32) 



For T <C Tc the quasiparticle density is suppressed by 
the supcrfluid pairing, so the scattering rate is reduced 
below this estimate. 

The results in figure [3] are accurate whenever F ^ cj 
and a; ^ A. In the BCS regime of weak coupling, with 
Tc ~ A <C /i; these conditions can both be satisfied for 
T ^ Tc where the estimate ([5^ applies. 

When a; ^ F, notably for the limit — !• of a static 
perturbation, the quasiparticle collisions serve to damp 
the countcrpropagating currents A J = ~ J\^. This 
causes the zero-frequency response to vanish for all tem- 
peratures. At small non-zero frequency, this damping will 
lead to a dissipative spin current, with A J cx iwAA/F. 

The fact that for a; — >■ the response is zero even for 
T > in the presence of quasiparticle scattering is very 
important if one wishes to measure the supcrfluid density 
in the way outlined in section |TT1 For that method to 
work one requires the response of A J to a perturbation 
A^ to be zero. Figure [3] indicates that in the absence 
of quasiparticle collisions x(w, T) becomes appreciable, a 
few percent of the maximal (non-interacting) value, for 
^ > |. On the other hand, as long as we can assume 
that F is much larger than uj, quasiparticle collisions will 
serve to ensure that A J = 0. Thus if a vector potential is 
applied to only one species, the spin current will be zero 
and the total current can be used to infer the superfluid 
density. 

This suppression of AJ at low temperature can be 
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understood as a form of spin drag. A vector potential 
A-f acting on the spin up particles causes a current 
which in turn causes a current Jj^, eventually (at low 
enough frequency) leading to a suppression of AJ. In 
non-superfluid ferromagnetic Fermi gases the spin drag 
has been investigated in Ref. [2^. The main difference 
in the superfluid systems we are considering here is that 
at w = there is no relaxation of the spin drag. This 
is because unless Cooper pairs are broken at w 7^ 0, the 
constituents of the Cooper pairs cannot move indepen- 
dently. 

V. EXPERIMENTAL CONSIDERATIONS 

The most natural experimental implementation of 
these ideas involves a modification of the geometry de- 
scribed in Ref. [T^ to the two-species Fermi gas. Thus, 
we consider the two species confined to a ring-like trap. 
The application of coherent o ptic al fields, using beams of 
non-zero angular momentum|15| . allows azimuthal vec- 
tor potentials to be imprinted on one, or both, of these 
species. Clearly the simplest case is to dress just one 
species (say spin-f), giving A-j- ^ and = 0. The 
methods described in Ref. J^, then allow the average 
azimuthal momentum of the dressed species, and hence, 
its current J^, to be measured spectroscopically. Tem- 
poral modulation of the amplitude of these optical fields 
and/or the detuning allows A^ to be time-dependent, al- 
lowing the use of this technique to probe the response at 
nonzero frequency. 

In order to understand the consequences of the results 
described above, we express the vector potentials in terms 
of the symmetric component A = 7^{A^ + A^) and a 
spin asymmetric component A A = {A^ — Ai). When 
just one species is dressed, A^ = A and A^ = 0, these 
components are both non-zero A = ^A and A A = A. 

Wc have shown that in the low frequency limit w <C F 
the response to AA is zero for any temperature. This 
implies that the system only shows a current response 
to A. This induces a linear response of the superfiuid, 
with a current density ([T]) set by the superfluid density 
Ps- The non-equilibrium state of the superfluid caus- 
ing this steady state flow remains (meta)stable provided 
the induced velocity remains below the superfluid crit- 
ical velocity. For optical dressing of a single species, 
A = iA and AA = A, the resulting current density 
of this species is = j + ^Aj = = -^A. 

Thus, a spectroscopic measurement of the total current 
of this species Jf = Vjf, where V is the total volume, 
allows a direct measurement of the superfiuid density ps- 
The superfluid density of a harmonically trapped Fermi 
gas has been measured using the collective modes [2?} . 
The method we propose here can be applied in a wider 
range of geometries (it does not rely on the harmonicity 
of the trap) . With local imaging, it would also be able to 
probe the local superfluid density in different parts of the 
atom cloud. Furthermore, the method allows interesting 



additional information to be obtained from the response 
at non-zero frequencies. 

For w > r the response to AA is non-zero and becomes 
appreciable for a; ~ A, as shown in Fig. [3l For w < 
2A but not oj A the response of the spin current 

— to AA provides a measure of the polarisability 
of individual Cooper pairs according to eq. p9)) . In the 
absence of a vector potential Cooper pairs are spherically 
symmetric. A time- varying vector potential AA acts like 
a time-varying electric field coupling differently on the 
two components of the Cooper pairs. This leads to a 
time-dependent polarisation of the Cooper pairs and thus 
to a time-dependent spin current. For ui > 2A the gauge 
field AA has enough energy to break Cooper pairs apart 
and the gas shows a response similar to that of a gas of 
non- interacting particles. 

VI. CONCLUSION 

We have investigated the response of a two-component 
Fermi gas to a vector potential which couples differently 
to the two spin species. The vector potential can be 
decomposed into a component acting in the same way 
on the two spin species (the spin-symmetric component) 
and a component acting in the opposite way on the two 
spin species (the spin-asymmetric component). We have 
shown that in the limit w — > the response to the spin- 
asymmetric component can be neglected, where uj char- 
acterises the time dependence of the vector potential. 
Thus only the response to the spin-symmetric compo- 
nent remains from which the superfluid density can be 
deduced, similar to what is discussed in Ref. [l^. We 
have also addressed the response of the spin current to 
the spin- asymmetric component for larger values of w. 
For UJ ^ 2A the response can be described in terms of 
a polarisability of the superfiuid, arising from the dis- 
placement of the average position of the spin up and 
spin down particles in the initially spherically symmetric 
Cooper pairs. This polarisability can be related to the 
strength with which the Cooper pairs in the Fermi gas 
arc bound. For cj ^ 2A the response of the system is the 
same as for a non-interacting Fermi gas, due to the fact 
that the perturbation has enough energy to overcome the 
Cooper pair binding energy. The response of the system 
is maximal when uj = 2A at which point the perturba- 
tion couples resonantly to the Cooper pair pair breaking 
process. Our results show that probing ultracold Fermi 
gases using vector potentials is a fruitful way forward in 
the study of fermion many body physics. 

Acknowledgments 

We acknowledge useful discussions with Z. Hadzibabic, 
N. Tammuz and E. Eliel on the experimental realisabil- 
ity. This work was supported by EPSRC Grant No. 
EP/F032773/1. 



[1] C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. 

Mod. Phys. 82, 1225 (2010). [17 
[2] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). [18 
[3] M. Zwierlein, A. Schirotzek, C. Schunck, and W. Ket- 

terle. Science 311, 492 (2006). [19 
[4] G. Partridge et al, Science 311, 503 (2006). 
[5] A. Leggett, Modem Trends in the Theory of Condensed 

Matter (Springer Verlag, Berhn, 1980). [20 
[6] K. M. O'Hara et al, Science 298, 2179 (2002). 
[7] J. GubizoUes et al., Phys. Rev. Lett. 91, 240401 (2003). [21 
[8] G. Regal, G. Ticknor, J. Bohn, and D. Jin, Nature 424, 

47 (2003). [22 
[9] S. Jochim et al, Phys. Rev. Lett. 91, 240402 (2003). 
[10] M. W. Zwierlein et al, Phys. Rev. Lett. 91, 250401 

(2003). [23 
[11] G. Ghin et al, Science 305, 1128 (2004). 
[12] J. Kinast et al., Phys. Rev. Lett. 92, 150402 (2004). [24 
[13] M. Zwierlein et al. Nature 435, 1047 (2005). 
[14] J. Dalibard, F. Gerbier, G. Juzeliiinas, and P. [25 
Ohberg, Artificial gauge potentials for neutral atoms, 
arXiv: 1008.5378, 2010. [26 
[15] N. R. Cooper and Z. Hadzibabic, Phys. Rev. Lett. 104, 

030401 (2010). [27 
[16] S. T. John, Z. Hadzibabic, and N. R. Cooper, Phys. Rev. 



A 83, 023610 (2011). 

S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. 
Pliys. 80, 1215 (2008). 

J. D. Jackson, Classical Electrodynamics, 2nd ed. (Jolm 
Wiley & Sons, New York, 1975). 

D. Forster, Hydrodynamic Fluctuations, Broken Symme- 
try, and Correlation Functions (W. A. Benjamin, Read- 
ing, Mass., 1975). 

A. Altland and B. Simons, Condensed Matter Field The- 
ory (Cambridge University Press, Cambridge, UK, 2006). 
P. D. Gennes, Superconductivity of metals and alloys (W. 
A. Benjamin, New York, 1966). 

L S. Gradshteyn and L M. Ryzhik, Table of Integrals, 
Series, and Products, 5th ed. (Academic Press, London, 
1994). 

A. Recati, 1. Carusotto, C. Lobo, and S. Stringari, Phys. 
Rev. Lett. 97, 190403 (2006). 

P. Nozieres, Theory of Interacting Fermi Systems 
(Addison- Wesley, Reading, MA, 1997). 
P. Nozieres and D. Pines, The theory of quantum liquids 
(Perseus Books, New York, 1966). 

R. A. Duine, M. Pohni, H. T. G. Stoof, and G. Vignale, 
Phys. Rev. Let t. 104, 220403 (2010). 
S. Riedl et aZ.. larXiv:0907.3814 Y2. 2009. 



